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I. INTRODUCTION 


The problems related to thin structural bodies (rods, beams, 
plates,and shells) with other bodies have widespread 
application in various fields of science and technology. The 
physical phenomena involved in the impact event include 
structural responses, contact effects and wave propagation. 
The problems associated with these are always topical 
issues in the field of applied mechanics. Since these 
problems belong to the problems related to dynamic contact 
interaction, their solution is connected with cumbersome 
mathematical tasks. To this end , several researchers had 
worked and some are still working on the dynamic behavior 
Ambartsumian 
examined the five fundamental differential equations 
describing the equilibrium of an orthotropic plate with a 
cylindrical anisotropy for the case when all radial planes 


of orthotropic rectangular plates. 
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Abstract— This work investigates the vibration of orthotropic 
rectangular plate resting on a variable elastic Pasternak foundation 
under the action of moving distributed masses. The governing equation 
is a fourth order partial differential equation with variable and 
singular co-efficients. The solutions to the problem are obtained by 
transforming the fourth order partial differential equation for the 
problem to a set of coupled second order ordinary differential 
equations using the technique of Shadnam et al[11] which are then 
simplified using modified asymptotic method of Struble. The closed 
form solution is analyzed, resonance conditions are obtained and the 
results are presented in plotted curves for both cases of moving 
distributed mass and moving distributed force. 


crossing the axis of anisotropy are the planes of elastic 
symmetry. Sveklo [2] suggested the contact theory for two 
anisotropic bodies under compression according to which 
the contact pressure is distributed over an elliptical contact 
region. The same structural effects are also true of the 
concrete slab in a composite girder bridge, but the steel 
orthotropic deck is considerably lighter, and therefore 
allows longer span bridges to be more efficiently designed. 
Awodola [3] studied the effect of plate parameters on the 
vibrations under moving masses of elastically supported 
plate resting on bi-parametric foundation with stiffness 
variation.Szekrenyes [4] investigated the interface fracture 
[1] in orthotropic composite plates using second order shear 
deformation theory. Yan [5] proposed elastic orthotropic 
models and used these in the nonlinear analysis of concrete 
structures subjected to monotonic or pseudo dynamic 
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loading. Since these models can appropriately describe the 
strain softening behavior of concrete beyond the peak stress 
and show good agreement with the strength envelope 
obtained from experimental results Hu and Yao [6] studied 
the vibration solutions of rectangular orthotropic plates by 
symplectic geometry method. In the same vien, Alshaya, 
Hunt and Rowlands [7] examined stresses and strains in 
thick perforated orthotropic plates. Gbadeyan and Dada [8] 
found the natural frequency of rectangular plates traversed 
by moving concentrated masses. Awodola and Adeoye [9] 
investigated the behavior of simply supported orthotropic 
rectangular plate by applying the technique of variable 
separable. Adeoye and Awodola [10] studied the dynamic 
behavior of orthotropic rectangular plate ^ with 
clamped-clamped boundary conditions by making use of 
the technique of nna Due to inability of researchers to 


Daga Wai + 2B gis 


W(x, y,t) + —— 


s ot? = ôt? 


12x + 3x?) -W (x, y, t) — S(12 — 13x + 6x? + DE 


W(x, yt) + Dy 2 z 
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solve orthotropic plates problems by analytical methods, 
this work aims at solving the governing equation by 
analytical solution and also considers the effect of the 
flexural rigidities in both x and y directions. 


II. GOVERNING EQUATION 

The dynamic transverse displacement W(x,y,t) of 
orthotropic rectangular plates when it is resting on a 
bi-parametric elastic foundation and traversed by 
distributed mass M, moving with constant velocity c; 
along a straight line parallel to the x-axis issuing from point 
y = s on the y-axis with flexural rigidities D, and D, is 
governed by the fourth order partial differential equation 
given as 


W(x,y,t) tuS 5 W, y,t) — phRo 


ot? 


W (x, y,t)] + Ko(4x — 3x? + uds y,t) + a 13 + 


"W(x, yt) +5 WO, y,t)] (2.1) 


SN, [M,gH(x — ct)H( — s) - M, (ÈW (x, yt) + 26 Way, D +c? 
Z W(x, y, t))H(x — crt)H (y — s)] = 0 


where D, and D, are the flexural rigidities of the plate along x and y axes respectively. 


3 

D Exh? D Eyh 

* 12(1-vxvy)' y 12(1-vxvy)' 
xVy xVy 





(2.2) 


E, and E, are the Young's moduli along x and y axes respectively, Go is the rigidity modulus, v, and vy are Poisson's ratios 
for the material such that E,vy = Eyv, , p is the mass density per unit volume of the plate, h is the plate thickness, t is the 
time, x and y are the spatial coordinates in x and y directions respectively, R, is the rotatory inertia correction factor, Kọ is the 
foundation constant, S, shear modulus and g is the acceleration due to gravity, H(.) is the Heaviside function. 


idi Dn (2.1), one obtains 


LÀ < y, t) + WZW (x, y, t) = phRo[—5 W (x,y, t) + ——W (x,y, t)] — 2B 


x: A 0y? Tx 


m ;WQo yt) - Dy a “wx, y,t)— D, = “we, y, t) + uw2W (x,y, t) — Ky (4x — 


3x? +x3)W(x,y, i + G,(—13 + 12x + 3x?) -W (x, y, t) -Gy(12— 13x + 6x? + x?) | (Q3) 


9? 
“wx, yt) tI EWG, y,t)] + Xr [M, gH (x — c,t)H (y — s) - M, Ga WQoy,t) 
26, Wo, ote ZWE y HG — eHG - s)] 
Simplifying ee (2.3) further, one obtains 
ðt 0^ 2B 
— pa WC y,t) + add y, pu = Ar t aa WO YO + Saga Wy, t)) — 
Do t UM 2 x = 
DD zW(x,y,t) — i s “wx, y,t) i tne yt) +w,W (x, y,t) as - 3x? 
+x3)W (x, e gem "m 13 + 12x + 3x?) Ž —W (x,y, 0-52 —13xt6x?tx es z (2.4) 
W (x, y,t) tI Swe A + Èr- Cs -gH(x — cyt)H(y — s) ur We, y,t) + 2c, 
Z way, 9 tc “wey, tH- crt)HY - )] 
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where (02 is the natural frequencies, n = 1,2,3,... 


The initial conditions, without any loss of generality, is taken as 


ô 
W(x, y,t) 0-7 xWQoy,t) (2.5) 


III. ANALYTICAL APPROXIMATE SOLUTION 


In order to solve equation (2.4), one applies technique of Shadnam et al which requires that the deflection of the plates be in 
series form as 


W x, y,t) = Xni Va Go y)Qs(t) (3.1) 


where 


P(x, y) = Vim (x) Pam) 
YPimX = Sin jmxX + AjmCOSGimxX + Bymsinh jx + CjmCOSHG jn X 


Pam = SiN@rmY + AnmCOSPnmY + Bamsinh@nmy + ChmcoshPnmy 
(o jm m m 
Cim = ua Qnm — r 


The right hand side of equation (2.4), taken into account equation (3.1), written in the form of series takes the form 


Dies RE V. 0 )0s 0 + 25 ENG- 3 DES r nEaN - 27. 
Wn (2, Ont) — BZ Va, XX) + o£. (x, Dt) -É (4x — 3x? x3) 9. (x,y) 
Qa (t) qs "m 13 t 12x + 3x?) = -Y n(X, y) Qn (t) — ^2 — 13x + 6x? + EÉ ay) (3.2) 
On(t) + EPC YOn) + XI. CE gli - 4080 -= s) - "E (9, (x, y)O. (0 + 
2c 2, (x, y)On(t) c? 22. Qn ()H c — eH O — )] = Xa V. (x, 3)8, (0 


Multiplying both sides of equation (3.2) by VU, (x, y) and integrating on area A of the plate and considering the orthogonality 
of 'P, (x, y), one obtains 


«(D = Bia f M ond RUN C ME ag ona) 


Qn (t) - 2 27 w, (x, y)Q, (t) - 2, w, (x, din + oV, x, "eo = 7o (4x = 


u mh u 20 
3x? + x3)W, (x, y)Q, (t) + 2 "i 13 + 12x + 3x?) Ž = Pn Go y)Q«(t) — 9 a2 — 13x 
(3.3) 
+6x° + xw, (x, E tI oY n(x, 399. E) FO (E n — c,t)H(y — s) 
-" C9 Gy), (t) 2c 2-9. (x, re +c? W(x, 9) Qn(t))H (x — c OHG — s) 
)] Pm (x, y)dA 
and zero when n zm 
where 
0* = f, P(x y)dA (3.4) 


Making use of equation (3.3) and taking into account equation (3.2), equation (2.4) can be written as 
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V, Go y) o0. (t + 0,(0] = BEM yz, [. [Ro A ws, Gy), (t) € 202 
n G0.) Aste Y E UE - 27309? w, (x, y)Q, (t) NU 
V. y)Q,(t — 7o (4x — 3x? + x3)W, (x, y) Ym "m y)Q,() + & "2-13 + 12x + 3x?) Sa) 
W(x, Qq(t) - -— — 13x + 6x? + x33) P 1n DP Ne 35:6 Es y y) $9) 


Qq(t)) + Xr (= - OF m(%, y)H (x — ¢t)H(y — s) - "E P, Qo y) V (HY) Gq CE + 2c, 


29409 ay Q5 y) 04) +c 29999 o, (x, Qq (0) H(x — c, H( — )))dA 


On further simplification of equation (3.5), one obtains 


On(t) + 9204 = SDR f, [Ro 


Dy a* sag 
H 


a? ay 0? ID 


W(x, y) 0.) + Maw, (x, yo (0) 
IP w, (x, y)Qs(t) — Vs (2, Y) Ct) - ME (x,y) 

Qq(t) + w2¥ (G9, 6:0,(0 — s - 3x?  x*)9, Qo y) Yn y) Qut) +È (-13 

a ea wx, Y) (E) 

V (VVU O) + XI. Cg (YH = c HQ — s) - = (Hy Qs y) Ws Qs y) 


ea”) e, (x,y) Ct) + TM Pm Y)Qq (OH E — DHO — s))]dA 


(3.6) 


Ti2x-R3x7)— — Yn (x, ef) — "(12 — 13x + 6x? + x?)( 


qa a? p y) 


T + 2c, 


The system of equations in equation (3.6) is a set of coupled ordinary differential equations 


where H(x — c,t) and H(y — s) are the Heaviside functions which are defined as 





HC = crt) = (for rece» HO- 8) = por yes Gm 
With the properties 
EE d 
(i) Ux -e0] 28(x—-et, gH) 2 66 — s) (3.8) 
T j 2crt ; 2 

(ii) f(x) (x — ct) 200 T 20, F(A — 8) = HOO (3.9) 

Using the Fourier series representation, the Heaviside functions take the form 
H(x —ct) o i4 iy, BORO O<x< (3.10) 
H(y-s)-i-iyk, 9800 9 Qe yc] 6.11) 





On putting equations (3.7) to (3.11) into equation (3.6) and simplifying, one obtains 
On(t) + Hos — Xa [ROS - i - 3 T5Qq(t) - 5T, 0.) + (wh = 
7 Fé)Qq(t) +Ë Tg + Tj) Q4 — Efa “E (Ts + Qoa EL SE Ly. Ez 
snes oye : E: MS _ ys. A E: piee ya TN -ER " E: “eee 


E sin(2j*1)rc,t x COS(2k-1)zts x sin(2k+1)7s 
ey te Ora eR DID 


mesa x Sin(2j*1)mc,.t x cos(2k+1)7s 
2c, (Ts + E35 ———————— Ej) ————— —— Ei ———— 
rC 9t m ZR -1 2j41 =) =1 2j+1 Èk- 1 2k+1 
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x sin(2k+1)7s x Cos(2j*1)nc,t x Sin(2j*1)mce,t 
— Yee Eiz: uu 5t = yea Fis a Tram AEn a? 


2k41 2j41 


* Por x sin(2k+1)7s 
+= (ER 1 Eis Bi pde om Deni Eig EE —————)»OQq(t) +c? Tro + 


cos(2j+1)mcrt x sin(2j+1)mc;t x cos(2k+1)ms x sin(2k+1)7s 
— un. Bg oe ed Big et E20 ~oa 





(3.12) 


2j41 2j+1 2k+1 

x COS(2j*1)nc,.t x Sin(2j*1)mc,t x cos(2k+1)7s 

)ttOR 1 Bag it Fx gu 2*1 Qua Eds ar a 

x sin(2k+1)ms co Mrg 
-YXj. Ej I) 0,0))] = X22 Dhar 228 Mp (OW G) 
which is the transformed equation governing the problem of an orthotropic rectangular plate resting on bi-parametric elastic 
foundation. 
where 


= f, [5% Go) Go y) + E Pa Go y) Vs Qo y)]d4 
T, = f, 2S SW Qo y) Go y)dA 
T, = f, 2 9 Go y) Is Qo y)dA 
T, = f, 2 [9 Go y) Qo y)dA 
Ff = f, V oy) (x, y)dA 


Ts = 4U, = 3U; + U3, Te = —13A, + 12A5 + 3A3 


T; = 12f, — 13f; + 6f3 + f4 + 12fs — 13f, + 6f, + fa 


T, = — f, VG, y) 9, Go y)dA 


by 
d 


t= f, Pa y)Us Go, y)sin(2j + 1)1xdA 


E = Í; V, (x, y)Ym (x, y)cos(2j + 1)nxdA 
E3 = f, V, Qc, y) Qo y)sin(2k + 1)mydA 
Ej = Í; V, (x, y)Ym (x, y)cos(2k + 1)rydA 


Eg = Ef, Eę= E}, E} =E}, Eš = Ey 


Ty = 2y (XY) Pm (x,y) dA 


16 lee Ox 
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(3.13) 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


(3.18) 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


(3.24) 


(3.25) 


(3.26) 
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Ej = Lx 5; (Hg Go Y) V, (x, y)sin(2j + 1)mxdA (3.27) 
Ej, = f, = (P, Qc y) Ym (x, y)cos(2j + 1)nxdA (3.28) 
Ej, = f, = P Q5 y) Vs Qo Y)sin(2k + 1)mydA (3.29) 
Ei, = if Zy, (x, y) V, (x, y)cos(2k + 1)rydA (3.30) 
Ej, EL Ei, = Et, Eis = Ej, Ete = Ef (3.31) 
Tio = 2- f, E OEY) En y)dA (3.32) 

Ei; = Í, Zo (x, yY))Ym(x, y)sin(2j + 1)txdA (3.33) 
Eis = f, DS P Qo y) Qc y)cos(2j + 1)xdA (3.34) 
Eio = Í, Pa y)Ym(x, y)sin(Zk + 1)rydA (3.35) 
Ežo = S, Zw, (x, y)) Wa (x, y)cos(2k + 1)nydA (3.36) 
Ej = Et, Ely Ed Ets = Ets, Ekg Ei (3.37) 


Wm (x,y) is assumed to be the products of functions Pym (x)Wp,, (y) which are the beam functions in the directions of x and y 


axes respectively. That is 


Wax y) = Von) Wn Cy) (3.38) 
where 
®,,(x) = sin? + Am cos™ + Bm sinh ^^ + C,,cosh —— P (3.39) 
®,,(y) = sin + Am cos L2 + Bn sinh %22 » + Cn cosh 27 (3.40) 
where Apm, Bpm Cog, Apm, Bom and "M are constants determined W the — conditions. And Vp, and V, are 


called the mode frequencies 
where 


Au 59m i= (3.41) 


Ly Ly 


Considering a unit mass, equation (3.12) can be re-written as 
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On (t) + @RQn(t) — 5: Laer [Ro Og (t) — = T204() RO - == TyQq(t) + 


* K x CO TC 
(ORFE —“2T5)Qq(t) + 9 (Te + Ty) Qq(t) — ap (Te + (Lf E; SO 


2j+1 
x sin(2j+1)mct x COS(2k-- 1)rs x sin(2k+1)7s 1 
-Eja Eu smeo ye 1 E3 NN pba— o 4T 
E: cos(2j * 1)zct Ei sin(2j+1)mct * cos(2k*1)ns — 
Èj- =1 Ë 2j+1 - Èj- =1 E ae ) + „k 1 Ez 2k+1 
x sin(2k+1)7s x COos(2j+1)rct P 
Èkr=1 Eg eO (t) + 2c(Ts + ai ZR Bil 2 - = $i Exo 
2k+1 2j+1 
sin(2j+1)mct x cos(2k+1)7s x Sin(2k-1)zs 
— anu OQ Eno au. 0 Àka En er) QE 1 Eis 


cos(2j-* 1)nct x sin(2j+1)mct x cos(2k+1)7s co * 
— iA Ea te Ej, ———>— E 
ou ye Dr dM UU — Yita Eis 
sin(2k+1)7s 


x COS(2j-*1)mct x sin(2j+1)mct 


2 = 
T ————))Qq(t) + c (Tio + (ap 1 Ei; 2j+1 X= 1 Eis 2j+1 
« cos(2k+1)ms x Sin(2k41)ms x cos(2j+1)mct 
ya, Ej, EU yp, py eS) + E re, pz, Ort 


sin(2j+1)mct x cos(2k+1)1s x sin(2k+1)zs 
Eh— 1 Qv E3, — —— — — E, — 
—Epjsi Zu. ?) = OR. 1 "ET Mia UE 





y E (3.42) 
)Qq(t))] = X2 nor Fm (CE) Pim (s) 
equation (3.42) is the fundamental equation of the problem. where 
M 
a= T p = LL, (3.43) 
V, (ct) = sinygt + Aj cosygt + Bmsinhymt + C,coshyst (3.44) 
Us) = sinv Ag cosv, + Bmsinhvm + C4coshvg (3.45) 
mE Pms 
= b V. = — 3.46 
Xm m PAS, an (3.46) 


3.1 Orthotropic Rectangular Plate Traversed by a Moving Force 


In moving force, we account for only the load being transferred to the structure. In this case, the inertia effect is negligible. 
Setting w = 0 in the fundamental equation (3.42), one obtains 


Ön (E) + (1 - Don CE) — A (URT, On (t) — 2BT2Qn(t) — DyTsQn(t) — DT Qs CE) 
-K,T;Q.(t) + ca (T; + T,)Qu(t)) — Y [RoT Gq (t) — 2BT2Q,(t)—DyTsQ,(2) — 40 
—DyFEQq(t) + (Mea Fe — KoTs)Qq(t) + Go(Ts  T7)04(] = Leer Elea M2 Yn (c£) Ps S) 

which can further be simplified as 


Qn(t) + On ni t (LRoT, Q, (t) B 2BT2Qn(t) ~ DyT3Qn(t) = D4T4Q4(t) U KoTs 
Qn Ct) + Go (Te + T7)Qn(t)) — *' Digna [URT Q4 (t) = ae () — D,T3Qq(t) — D;T, (3.48) 
Qq(t) + (uo, Ft — KoTs)Qq(t) + Go(Ts  T7)Q4(0] = Xg2i X21 Mgt ^W, (ct), s) 


On expanding, re-arranging and simplifying equation (4.48), one obtains 


Önt) +É Qu t) + ;Eseuaen (RoT G4 Ct) — 2BT2Q4 (t) — Dy TaQa (Ct) 


[1-T*uRoT1] [1-T*uRoT1] 
T*Mg 


—D,T4Qq(t) (uo Fz — KoTs)Qq(t) + Go(Ts + T7)Qq(t)) = uR F] 


(3.49) 
Pm (ct) "Pm (S) 


where 


www.ijaers.com Page | 32 


Awodola T.O et al. International Journal of Advanced Engineering Research and Science, 8(6)-2021 


* 1 Fy * 
T = Tg” Xi-ü- 900. Jé = —2BT, — DyT3 — D,T, — KoTs + Go(T> + T7) (3.50) 


: ; : Y ; 
For any arbitrary ratio Y, defined as t* = iy one obtains 


T* 





Y=- eq o(t*?)+... (3.51) 
For only o(t*), one obtains 
T'-Y (3.52) 
Applying binomial expansion, 
1 = 14 YuRSF; + o(Y2)4... (3.53) 





1-YLRgF, 
On putting equation (4.53) into equation (4.49), one obtains 
Qn(t) + Q — Y) + YURT, + o(Y?)4...)0,(t) + Y(1 + T*URoT, + o(Y2)4...) 
Yi-iaen (URoF{Zq(t) — 2BF; Q4 (t) — DyF3Qq(t) — D;FzQq(t) + (u$ — Ko)F§Qq(t) + (3.54) 
Go (F5 + F7)Q4(t)) = (1 + YuRoT, + o(Y¥7)+...)Mg®,,, (ct) 9, (s) 
Retaining only o(Y), equation (4.54) becomes 
Ön CO + Q1. YHRoT1) — YJ6)Qn(t) + Y Lgeaqen (4RoF1Q4(t) — 2BF2Qq(t) — DF; 


A . : . (3.55) 
Qq(t) — Dy Fi Qq(t) + (ue — Ko)F2 Qa CE) + Go(Fe + F7)Qq(t)) = YMg WW, (ct) Ws) 
which is simplified further as 
Qn(t) t J7Qn (©) F Y cin (URF? Q,(t) = 2BF;Q, (t) NE Dy F3Q,(t) 7 D,F4Qq (t) T (3.56) 
(Hag — Ko)Fs Qq(t) + Go(Fé + F7)Qq(t)) = YMg Pm (ct), (s) 
where 
J; = x& (d + YuRoT;) — We (3.57) 
Using Struble's technique, the solution to the homogeneous part of part of equation (3.57) is assumed to take the form 
Qn(t) = q(n,t)cos(yst — p(n,t))+... (3.58) 
where 
p(n, t) =e, (3.59) 
and 
xñ- 
$(n,t)- m + ly (3.60) 
On putting equations (3.59) and (3.60) into equation (3.58), one obtains 
2_7* 
Qn (t) = &cosQy,t — GERE 1) (3.61) 
On further simplification, one obtains 
Qa (t) = £ycos(v,t — tn) (3.62) 
where 
2. 7* 
Un = Xn - CEE) (3.63) 


is the modified frequency for moving force problem for orthotropic rectangular plate resting on variable elastic bi-parametric 
foundation. 


Using equation (3.63), the homogeneous part of equation (3.56) 
On(t) + vi Q4 (t) = 0 (3.64) 


Hence, the entire equation (3.56) gives 


On(t) + v2Q, (t) = YMg¥ nm (ct) m (s) (3.65) 


Re-writing equation (3.65), one obtains 
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On(t) + v2Q,(t) = YMgW,, (s)[siny; (t) + AmCOSXmt + Bmsinhymt 


3.66 
+C,,coshy,,t] (3.66) 


To obtain the solution to equation (3.66), one makes use of Laplace transformation techniques to obtain 
= en > 


(uie — a — Bn (x, — v2) (Qmsinuy,t — v_zsinhy mt) + C, v, (x2, — v2) (3.67) 
(coshy,,t — cosu,t)] 


which on inversion yields 
MgYPm o 


W(x,y, t) EX m=i mci v n(X&-v ja zi + Un) Gn sinv,t = UnSiNXmt) ~ AmUn (Xm 
+u2)(cos%mt — cosu,t) — "pd — v2)(%msinu,t — v nsinhtmt) + Cav 2, 


—v2)(coshy,,t — cosu nt)](sin = bmeg Ajmcos P x + Bjmsinh 7 Pimi pues) 


$ jm 


jm E 











(sin $m "— + AnmCcos 777 Pum , + -— d + Chmcosh fam yy 


which is the transverse dem response to a moving force of orthotropic rectangular plate resting on variable elastic 
bi-parametric foundation. 


3.2 Orthotropic Rectangular Plate Traversed by a Moving Mass 


In moving mass problem, the moving load is assumed rigid, and the weight and as well as inertia forces are transferred to the 
moving load. That is the inertia effect is not negligible. Thus w + 0 and so it is required to solve the entire equation (3.42). To 
solve the equation,one employs analytical approximate method. This method is known as an approximate analytical method of 
Struble. The homogeneous part of equation (3.42) shall be replaced by a free system operator defined by the modified 
frequency vn. Thus, the entire equation becomes 


A * N00 * * 2j+1 t oo * 
On(t) + v2Qn(t) + ap" Xi CF + OR E; EDD ye E; 


2j+1 

cm OQ. , f; EEUU c yt, Et EQ 
mom yy pebioi op, ge eU" xe n 
RG g(t)  2c(F3 + (Eu. ES SECO ym, p; E 
YQ, ri SEO ya, pi OS) ue, gi tope 


Ej, sin(2j+1)rct x cos(2k+1)7s x Sin(2k+1)rs 
- tI pj, REUS L p M 
ij E 2j41 eg ev 2k+1 Ya 2k+1 


* 2j+1)nct x 2j4*1)nct 
DÅ) + c (3, + Qr Ej, SSA L 2 ae 


2j+1 2j+1 
cos(2k+1)ms x sin(2k+1)7s x cos(2j+1)mct 
Et, = gi Smee pi EOS MESE 3.69 
Qs 2k41 Yi 2k41 ) + nT =1 2j+1 ( ) 
Eż, sin(2j+1)rct â pare rere x sin(2k+1)ms 
- tI pj, prm 
ij E 2j41 ta 9 ev 2k+1 Xia 2k41 





DUOI = X221 5. Pm(Ct) Pm (s) 


where p* = E 


On expanding,simplifying and rearranging equation (4.88), one obtains 


^ eas x 2j+1)rct T 2j41)nct 
On(t) + 2cap'(F += (Xj, E; SEPT Lye, pu PEO yea Et 


a 2j+1 
sae ye, eM e Log, rn AE Te, E, 
OLD + (Y, Fis esene yp 1 Eje 0727-0, (1) + 
(RC. — ap +a Ea BS - yg, E 
Or, nm EU yp, gene 4 ye, gita. 
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Ei sin(2j+1)mct x coS(2k+1)mct x sin(2k+1)1ct 
um 2. POR m oles a 


c^ap' (Fio + ZER 1 Er; n Pe 1 Eis an U 1 Eio 
— YQ p RO) + SR Ej ma - ym Er 
MOADE (iy pj TOU o yp, p; Rn (Ct) 
ap" Tena là + ay a. Ej SO ym, By (Di ES 
De — Yga Ej )tuQg nO yr E 
EG (Xp. p; E ya Bg y, (t) + 


cos(2j * 1)zct x Sin(2j*1)mct x Cos(2k+1)rct 
2c(F3 Eg——— Eio Ek — 
( 9 +40 =1 m gx T 2j41 eu Jus 1 2k41 


cos(2k+1)mct 
2k+1 


cos(2k+1)mct x sin(2k+1)1ct 


2Kk+1 


sin(2k+1)mct x COS(2j-*1)nct Ek. sin(2j+1)mct 

zZ Et,————— Ej}, — — — um s nr dca 
yu 1 2k41 )+ lcm 1 2j41 bs ABE 2j41 ) 

x cos(2k+1)1ct x sin(2k+1)mct 


2 * 1 oo * 
ae kii ))Q,(t) T c^ (Fio + za Qa Ej; 


« Sin( + cos(2k+1)mct co " 
m 1 Eig Eio I zki — = et E59 


TO Eis = 1 Eie 


2j*1)nct 
meee 


2j*1 
sin(2k+1)rct 


x cos(2j+1)mct « Ssin(2j+1)rct 1 
DTN (XR. Ej, SO L ym) Ej eet 


2j+1 i 4T 

(3.70) 
" 2k41)nct y 2k+1)rct 
EL Ej coset ye 1 Ej, 585227750, (t)] = Sn (ct) CS) 


2k+1 2k+1 


Applying modified asymptotic method of Struble, the solution to equation (3.70) takes the form 


where 


and 


Qa (t) = Y(n, t)cos(u,t — G (n, t)) + v4Q,(t)... (3.71) 
Y(n, t) = Bee Ist (3.72) 


1 IT x COS(2k-- 1)rts E: sin(2k41)rs 
F(n,t) = Ti (F5 + l0 1 Eaa ra eo 2) 





(3.73) 
bryk x Cos(2k+1)rct x Sin(2k+1)rct 
—c*^mq' (Fh + ONA 1 Ez 3 oe Èk=1 E2 — le ts 
On putting equations (3.72) and (3.73) into din (3.71), one obtains 
-ap* 1 nns 4 cos(2k+1 re 
Qn (t) = Ze"? ht cos[unt — (vap (Fe + E (Diver Ej SEO yp, Bg 
sin(2k+1)7s 2 &YIK x Cos(2k+1)rct _ x sin(2k+1)mct 3.74 
a) — ap y+ = (Sh Ej SO ey Ez DT) S) 
»t - Tr] 
On further simplifications, one obtains 
Qa (t) = Ze~%? Istcos[B,t — En] (3.75) 
where 
1 prrs 4 Cos(2k+1 pe Sinak+1 
B, Up Nd (+208, SS Cy, pg Em) 
n (3.76) 
eer x Cos(2k+1)rct x Sin(2k+1)rct 
—c*ap* (Fio + ar 1 E T T Eka E ey) 
is the modified frequency representing the frequency of the free system. 
Using the same argument used to solve equation (3.56), the homogeneous part of equation (3.70) becomes 
Qn(t) + PaRa (t) = 0 (3.77) 
Hence, the entire equation becomes 
On(t) + BaQn(t) = 2 On (ct), (s) (3.78) 
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Rewriting equation (3.78), one obtains 
QC) + B2Q,(t) = oo o, (s)[siny,, t + Amcosymt + Bmsinhymt + C,,coshy,t] (3.79) 


Following the procedures applied to solve equation (3.66), one obtains 


Qn(t) = a (GÀ, T B2) s sinf, t = PnSiNXmt) m AmP Wn + B) 
(cosy,,t — COSByt) — Bm (x2, — B2) (ra sinfst — Bnsinhymt) + Cs. Qr, — B2) (3.80) 
(coshy,,t — cosB,,t)] 


which on inversion yields 


oo oo Om H . 
Way, = Dimer Litmat gp cot op [Oth + Ba) msinBat — Basin t) — Anf (th + Ba) 
(COSXmt i cosp,t) C Ha z Br) (Xmsinpnt im iS n + Cmbn Xh — Be 


(coshynt — cosByt)](sin $E x + A5, cos x + Bj, sinh x + C, cosh x) (3.81) 


(sin Es Pn F + Ang cos =" Pici zu + By, sinh =" fin, " + Chmcosh fin y 
which is the transverse Menem response to a moving mass of an orthotropic rectangular plate resting on variable elastic 
bi-parametric foundation. 
3.3 ILLUSTRATIVE EXAMPLES 
3.3.1 Orthotropic Rectangular Plate Clamped at All Edges 


For an orthotropic plate clamped at all its edges, the boundary conditions are given by 


W(0,y,t) 20, W(L,y,t) 20 (3.82) 
W(x,0,t) 20, W(x,L ys t) -0 (3.83) 

OW(Oyt) _ OW(Lyy,t) _ 
mc m 0, t wem 0 (3.84) 

OW(x,0,t) _ OWQoLy,t) _ 
"ue cet x c9 (3.85) 

Thus,for the normal modes 

Spm(0) = 0, £5 (L4) = 0 (3.86) 
fgm(0) = 0, Som(Ly) = 0 (3.87) 


O85m(0) =0 O85m(Lx) 


am i às = 0 (3.88) 
O£qm(0) O£qm.y) m 
x m 0, x 0 (3.89) 
For simplicity, our initial conditions are of the form 
W(x,y,0) = 0 = War (3.90) 


ot 


Using the boundary conditions in equations (3.86) to (3.89) and the initial conditions given by equation (3.90), it can be shown 
that 


Apm m sinhéyg -sinépm 2 cosgnm—coshinm (3.91) 


cosgpm—coshinm singpm+sinhépm 
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A _ sinhggm—singgm | COSÉqm-coShéqm 
qm ^^ costo -coshíqm  sin€gmtsinhégm 





In the same vein, we have 


Jum sinhém—singm _ cosém—coshém 
m cosém—coshém sin£g +sinhëm 


Bom = —1, Bam 2-1, >Bm=-1 
Com = —Apm Cam = —Agm, > Cm 7 —Am 
and from equation (3.93) , one obtains 
cosé,,coshé,, = 1 


which is termed the frequency equation for the dynamical problem, such that 


č = 4.73004, č, = 7.85320, & = 10.9951 


(3.92) 


(3.93) 


(3.94) 


(3.95) 


(3.96) 


(3.97) 


On using equations (3.91), (3.92), (3.93),(3.94),(3.95) and (3.97) in equations (3.68) and (3.81), one obtains the displacement 
response to a moving force and a moving mass of clamped orthotropic rectangular plate resting on bi-parametric condition 


respectively. 


3.3.2 Graphs of the Clamped-clamped End Conditions 


Figures 1 and 2 display the effect of foundation modulus K, on the deflection profile of clamped-clamped orthotropic 


rectangular plate under the action of load moving at constant velocity in both cases of moving distributed forces and moving 


distributed masses respectively. The graphs show that the response amplitude decreases as the value of foundation modulus K, 


increases. 
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Fig.1: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying K, and Traversed by Moving 


Force 
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Fig.2: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Ky and Traversed by Moving 
Mass 


Figures 3 and 4 display the effect of shear modulus G, on the deflection profile of clamped-clamped orthotropic rectangular 
plate under the action of load moving at constant velocity in both cases of moving distributed forces and moving distributed 
masses respectively. The graphs show that the response amplitude decreases as the value of shear modulus G, increases. 
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Fig.3: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Gy and Traversed by Moving 
Force 
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Fig.4: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Gy and Traversed by Moving 


Mass 


Figures 5 and 6 display the effect of flexural rigidity of the plate along x-axis D, on the deflection profile of 
Clamped-clamped orthotropic rectangular plate under the action of load moving at constant velocity in both cases of moving 
distributed forces and moving distributed masses respectively. The graphs show that the response amplitude decreases as the 


value of flexural rigidity D, increases. 
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Fig.5: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying D, and Traversed by Moving 
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Fig.6: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying D, and Traversed by Moving 
Mass 


Figures 7 and 8 display the effect of flexural rigidity of the plate along y-axis Dj, on the deflection profile of Clamped-clamped 
orthotropic rectangular plate under the action of load moving at constant velocity in both cases of moving distributed forces and 
moving distributed masses respectively. The graphs show that the response amplitude decreases as the value of flexural rigidity 


D, increases. 
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Fig.7: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Dy and Traversed by Moving 
Force 
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Fig.8: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Dy and Traversed by Moving 
Mass 


Figures 9 and 10 display the effect of rotatory inertia R, on the deflection profile of Clamped-clamped orthotropic rectangular 
plate under the action of load moving at constant velocity in both cases of moving distributed forces and moving distributed 
masses respectively. The graphs show that the response amplitude decreases as the value of rotatory inertia R, increases. 
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Fig.9: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying R, and Traversed by Moving 
Force 
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Fig.10: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying R, and Traversed by Moving 


Mass 


Figure 11 displays the comparison between moving force and moving mass for fixed values of Ry, Go, Ko, D, and Dy. 
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Fig.11: Displacement Profile of Comparison between Moving Force and Moving Mass 


IV. CONCLUSION 


In this research work, the problem of vibration of 
orthotropic rectangular plate under the action of moving 
masses and resting on a variable elastic Pasternak 
foundation with clamped end conditions has been studied. 
The closed form solutions of the fourth order partial 
differential equations with variable and singular 
coefficients governing the orthotropic rectangular plates is 
obtained for both cases of moving force and moving mass 
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using a solution technique that is based on the separation of 
variables which was used to remove the singularity in the 
governing fourth order partial differential equation and 
thereby reducing it to a sequence of coupled second order 
differential equations. The modified asymptotic method of 
Struble and Laplace transformation techniques are then 
to obtain the analytical to the 


employed solution 


two-dimensional dynamical problem. 


The solutions are then analyzed. The analyses show that, for 


Page | 42 


Awodola T.O et al. 


the same natural frequency and the critical speed, the 
moving mass problem is smaller than that of the moving 
force problem. Resonance is reached earlier in the moving 
mass system than in the moving force problem. That is to 
say the moving force solution is not an upper bound for the 
accurate solution of the moving mass problem. 


The results in plotted curves show that as foundation 
modulus K, and the shear modulus G, increase, the 
amplitudes of plates decrease for both cases of moving 
force and moving mass problems. As the rotatory inertia 
correction factor R, increases, the amplitudes of plates 
decrease for both cases of moving force and moving mass 
problems. As the flexural rigidities along both the x-axis 
D, and y-axis Dy increase, the amplitudes of plates 


decrease for both cases of moving force and moving mass 
problems. 


It can be shown further from the results that for fixed values 
of foundation modulus and shear modulus, rotatory inertia 
correction factor, flexural rigidities along both x-axis and 
y-axis,the amplitude for the moving mass problem is 
greater than that of the moving force problem which implies 
that resonance is reached earlier in moving mass problem 
than in moving force problem of simply supported 
orthotropic rectangular plates resting on bi-parametric 
foundation. 
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